c program for curve fitting